function innerPoints = triangle_inner_points(V, n)
if n < 2
    error("n >= 2 required.");
end
[I, J] = meshgrid(1:n-1, 1:n-1);
valid = (I+J)<=n-1;
I = I(valid);
J = J(valid);
K = n-I-J;
mod = [K, J, I]/n;
innerPoints = mod*V;
end